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Abstract 

We reexamine the possibility of reconstructing the initial fluxes of supernova neutrinos emit- 
ted in a future core-collapse galactic supernova explosion and detected in a Megaton-sized water 
Cherenkov detector. A novel key element in our method is the inclusion, in addition to the to- 
tal and the average energies of each neutrino species, of a "pinching" parameter characterizing the 
width of the distribution as a fit parameter. We uncover in this case a continuous degeneracy in the 
reconstructed parameters of supernova neutrino fluxes at the neutrinosphere. We analyze in detail 
the features of this degeneracy and show how it occurs irrespective of the parametrization used for 
the distribution function. Given that this degeneracy is real we briefly comment on possible steps 
towards resolving it, which necessarily requires going beyond the setting presented here. 
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I. INTRODUCTION 



The detection of neutrinos from the SN1987A by the Kamiokande [l|, 0] and 1MB 
collaborations not only confirmed the basic picture of the core collapse supernovae (SN) but 
also initiated a whole new field of observational supernova astrophysics. In fact, the obser- 
vation of SN neutrinos is the unique way, except possibly for gravitational wave detection, to 
directly probe into the interior of a star which is undergoing gravitational collapse. Despite 
this pioneering observation, a precise understanding of the physics of a SN explosion still 
eludes us. Spherically symmetric models of iron core collapse do not explode, even with 

In the two-dimensional models 




111 ], reflecting the increased 



solid neutrino transport 0, 0, f| and general relativity 
the outcomes vary qualitatively and quantitatively 
complexity of the physics involved. 

Owing to their feeble interactions, neutrinos are able to escape from deep inside the star. 
Therefore SN neutrinos can provide us with information about the highly dense inner layers, 
where the SN explosion is initiated. Moreover, the composition of the SN core is such that the 
reactions involved in the creation and annihilation of neutrinos are different for the different 
neutrino flavors. These differences are imprinted in the emerging neutrino fluxes. These are 
non-thermal and can be characterized by the total energy emitted, the mean energy and 
the so-called pinching parameter which controls the width of the distribution. Therefore, 
a flavor-dependent reconstruction of SN neutrino flux is a useful probe for diagnosing SN 
core. (For early references, see e.g., 12|.) The physics which can be probed would include 
that governing matter under extreme conditions, such as information about the equation of 
state, and the explosion mechanism itself l^ . 

The potential of neutrinos in probing the SN core results, on the one hand, from the 
huge flux of them that will be released, corresponding to around 99 % of the total energy 
emitted in the SN explosion. On the other hand it is helped (though it may sound ironic) by 
the flavor mixing properties of neutrinos, including the large "solar" neutrino mixing angle 
obtained by KamLAND and the solar neutrino observation, in particular, SNO 15] . 
For an updated global analysis of the data of the various neutrino oscillation experiments, 
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see e.g., 

If a supernova explosion takes place in our galaxy the number of neutrino events expected 
in the current and planned neutrino detectors will be enormous 18|, ll9j]. Among all neutrino 
detectors water Cherenkov detectors are most likely the ones which can run long enough 
to watch galactic SN over long enough time scales. It is therefore important to establish a 
strategy for diagnosing the core of SN by using neutrino observation by water Cherenkov 
detectors. 
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Inverse beta decay u e +p — > e + +n provides the main neutrino detection reaction in such 
a detector. According to the standard picture of supernova neutrino propagation, except for 
the case of inverted mass hierarchy and large 613 (sin 2 # 13 > 10~ 3 ), the energy spectrum of 
i> e at the Earth is expected to be a strong mixture of the original v e and v x (we collectively 
denote z/ M , i/ T , u^, v T as v x since at first approximation their properties are identical and 
can be treated as a single species), modulated by the large value of solar mixing angle, #i 2 , 
well determined by solar and KamLAND data 16l Il7|. By performing a suitable simulation 

2G, L^lf how one can determine the 



of the high-statistics SN data it was shown in Refs. 
parameters of the original neutrino spectra, v e and v x in terms of the v e signal detected in 
inverse beta decay. However, this result has been obtained under the assumption that we 
know the parameter which describes the "pinching" of the SN neutrino fluxes. (See Sec. Ill Al 
for the pinching.) 



In this paper, we extend our previous work in 20j by including the pinching of the 
SN neutrino spectra as fit parameters 1 . In the extended framework, regrettably enough, 
we face with an important difficulty in reconstructing the original neutrino spectra. Most 
significantly, we find that there exists a persistent continuous degeneracy in the flavor- 
dependent determination of the luminosity and the spectra of SN neutrinos. 

In order to make our point and to indicate how serious the degeneracy problem is, we 
show in Fig. [1] examples of two degenerate v e energy spectra at the Earth (upper panels) as 
well as their difference divided by the true spectrum or fractional difference (lower panels) 
for two different parametrizations of the SN neutrino flux. It is remarkable that, despite the 
large difference in the primary neutrino spectra at the SN core particularly for u x , the true 
and fake spectra agree with each other within better than 1% level over a wide energy range. 
It is evident that the degeneracy is present already in the bare neutrino fluxes reaching the 
detector (before taking into account experimental uncertainties such as energy resolution) 
and it looks so perfect that it is unlikely to be resolved even with the extremely high statistics 
expected in a Megaton detector 0, 25, 26]. Hopefully, our negative result will stimulate 
further studies towards a full diagnosis of the SN core by neutrino observations. These 
include a better theoretical understanding of the neutrino fluxes formation, as well as an 
optimization of the information provided by the different reactions in the neutrino detectors. 

The plan of the paper is as follows. In Sec. [Ill we describe the basic features of supernova 
neutrinos, including the different parametrizations of the initial neutrino spectra usually 
considered in the literature, as well as the effect of neutrino flavor conversion before reaching 



1 A study of the relevance of this parameter in the analysis of the neutrino signal from SN1987A can be 



found for instance in Refs. 22, 12 
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FIG. 1: Examples of degenerate determination of SN astrophysical parameters. The two degen- 
erate neutrino energy spectra at Earth (upper panels) as well as their fractional difference (lower 
panels) resulting from neutrino fluxes with two different energy spectra at the SN core and the 
indicated parameter sets. The left and right panels correspond to the Garching parametrization 
and the modified Fermi-Dirac parametrization, respectively. Input values are the same ones shown 
in Table HI While it is not clear from the legend in the right panel, we note that the both true and 
fake values of average energies for the Fermi-Dirac parametrization give the same values as in the 
case of Garching one. Normal neutrino mass hierarchy has been assumed. Notations of neutrino 
species and the definitions of the parameters involved are given in Sec. HH 

the detector. In Sec. HHI we will specify the assumptions used in our analysis. Those include 
the particular values characterizing the initial neutrino fluxes, the details of the detector, 
and the neutrino mass scheme. In Sec. IIVI we demonstrate the existence of a continuous 
degeneracy in the determination of the astrophysical parameters. We discuss the robustness 
of this finding, which results neither from the particular parametrization taken, nor from a 
particular set of SN parameters. Finally, in Sec. [V] we speculate possible ways that could 
help overcoming the degeneracy problem. 
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II. SUPERNOVA NEUTRINOS 



A. Basic features of neutrino spectra from supernova 

To a crude approximation the proto- neutron star is a black-body source for neutrinos 
of all flavors. For the case of v e and v e the dominant reactions are the charged-current 
(CC) interactions with nucleons e~ + p <-> v e + n and e + + n <-> v e + p. The other flavors, 
z^, i/ r , v T which in this paper we shall collectively denote by u x , interact with the 
surrounding matter via neutral current (NC) interactions, e.g. Bremsstrahlung, neutrino- 
pair annihilation or neutrino-nucleon scattering. These processes keep u e , v e and v x in local 
thermal equilibrium up to the radii where these reactions become inefficient (neutrinosphere). 
Beyond these radii neutrinos freely stream. Taking into account the hierarchy in the cross 
sections, acc > a NC as weu as the richer neutron composition than protons, one expects 
the average neutrinosphere radii of v x , v e , and v e to obey r x < r s < r e , so that u x (u e ) 
decouples at the highest (lowest) temperature. This translates to an ordering of the average 
SN neutrinos (E Ue ) < (E D ) < (E„J), the exact degree of difference still under 



energies o 
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29]. See, for example, 30|, l31| for more about physics in the proto-neutron 



debate 
star. 

The location of the neutrinospheres does not only depend on the neutrino flavor but also 
on its energy. The energy dependence of the cross sections of the processes involved makes 
neutrinos with different energies decouple from the proto-neutron star at different radii and 
therefore different local temperatures. For this reason the spectrum of the neutrinos leaving 
the star does not present a thermal distribution. The possibility of reconstructing the flux 
parameters of the three effective flavors from observation would lead to a "neutrino imaging" 
of the proto-neutron star. 

There are different ways to characterize the non-thermal spectra of the neutrino fluxes 
emerging from the SN. Among them there are two parametrizations that have been exten- 
sively used in the literature. One is the Fermi-Dirac distribution motivated by the equilib- 



rium distribution of neutrinos inside the star 22j 



where E is the neutrino energy, and T Ua and r] Ua denote an effective temperature and degen- 
eracy parameter (chemical potential), respectively. The distribution is normalized so that 
stands for the total number of v a emitted. The function f n (r] Uci ) is defined as 
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The mean energy and the total energy released are consequently (E Ua ) = 
[h(Vu a ) / f2{Vu a )}T Ua and E^ = ® Va (E Ua ), respectively. Note that it must be understood 
that <& Vx does not refer to the sum of the flux of the non-electron species (despite that we 
treat them as a single species) but the individual one as follows, 

^ = $^ = $p m =^ t = $p t , (3) 

and so as for E^, throughout this paper. 

A way to determine how much a spectrum deviates from being thermal is to use the 



pinching parameter 27| defined as the ratio of the first two moments 



P-^i- (4) 



(Ey 

A spectrum that is thermal up to its second moment has p = PFD,r?=o ~ 1.3029, while p < 
Pfd,?7=o implies a pinched spectrum (high- and low-energy parts of the spectrum relatively 
suppressed) and p > Pft>,ti=o is an anti-pinched spectrum (high- and low-energy parts of 
the spectrum enhanced). For the Fermi-Dirac distribution with an arbitrary t] the pinching 
parameter p is related to t] as 

P = [U(v)f2(v)/fKv)} 

. In the left panel of Fig. [2] we show in solid red lines the explicit dependence of the pinching 
parameter on the effective degeneracy parameter t]. 

The curve becomes flat at negative rj, which reflects the fact that the function tends to 
the Maxwell-Boltzmann spectrum at rj — > — oo, which does not differ much from the Fermi- 
Dirac one, rj = 0. In dashed lines we present the strong dependence of (E) /T on 77 for 
pinched distributions. 

A second parametrization using the form of power times exponential has been recently 
introduced by the Garching group as a better parametrization of their simulation results 
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30 



In this case the parameters characterizing the distribution function are the total number 
of u a emitted, $^ Q , the mean energy, (E Ua ), and a parameter j3 Va which is related to the 
pinching parameter via 

p = (ft* + 1)/& . 

In the right panel of Fig. [2] we show in solid red lines the dependence of p on (3 Ua . In contrast 
to the Fermi-Dirac parametrization there is no asymptotic limit for the width (E 2 ) / (E) 2 , 
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FIG. 2: The Left panel gives the pinching parameter p = (-E 2 ) / (E) 2 (left label) and the ratio 
(E) /T (right label) for the Fermi-Dirac distribution represented as a function of r/ by the solid 
and dashed curves, respectively. The critical value of p which separates between a pinched and 
an anti-pinched distribution is indicated by the dotted line. The right panel gives the pinching 
parameter for the Garching parametrization as a function of (3. 



hence this parametrization can reproduce better anti-pinched distributions. In Fig. [3] we 
illustrate the effect of the pinching on three spectra with different /5's. One can see how the 
pinched spectrum is suppressed at low and high energies with respect to the non-pinched 
one. This in turn shows a similar suppression in comparison with the anti-pinched spectrum. 
We observe in the same figure how this behavior holds also for the case of the Fermi-Dirac 
distributions. The results shown above agree with those found in 30] where the similarities 
as well as differences of these 2 parametrizations were explored. 



B. Flavor conversion of supernova neutrinos 



The dynamics of neutrino flavor conversion in a typical iron-core SN can be factorized 
into two different parts 2 : the propagation through the inner layers, where the high neutrino 



2 There are also particular scenarios with very shallow electron density profiles, like the early stages of 
O-Ne-Mg core SNe 33[, where the dynamics in the inner and outer regions can not be decoupled. In 
this case collective effects and MSW resonances are not clearly separated, what can lead to interference 
effects both in normal and inverted mass hierarchy 34[. In the following, though, we will not consider 
this particular case. 
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FIG. 3: The solid curves show examples of anti-pinched {(3 = 2, shown in blue), un-pinched 
(j3 = 3.31, in black), and pinched (/? = 5, in red) Garching distributions normalized to unity. The 
mean energy has been set to 15 MeV. The corresponding Fermi-Dirac distributions for the same 
value of (E) and p are also shown by the dashed lines. Note that there is essentially no anti-pinched 
Fermi-Dirac distribution for (3 = 2 because of the saturation property of p shown in Fig. [2j 



density can lead to non-linear collective effects [34], [35|, [36|, [37|, [38|, [39|, |40|, |41|, |42|, |43|, |44( , 
and the evolution in the outer layers. Right above the neutrinosphere the neutrino density 
is so high that neutrino-neutrino interactions must be taken into account. The presence 
of neutrino self-interactions potentially give rise to collective phenomena on the neutrino 
propagation. These include synchronization, pair conversions of the kind v e v e —> v x v x during 
the so-called bipolar regime, and spectral split. These effects, though, arise only in the case 
of inverted neutrino mass hierarchy. 

In the outer layers the neutrino-neutrino interactions can be neglected and therefore it is 
enough to consider the Mikheev-Smirnov-Wolfenstein (MSW) effects induced on neutrinos 

3 Si. The neutrino flavor evolution can be then described in terms 



by background matter 

of two independent two-level crossings associated to the atmospheric and solar mass squared 



splittings 47|, |48j . Let us denote the former (latter) the H (L) level crossing. The relation 
is simplified by the fact that the L level crossing is adiabatic, given the confirmation of the 
large mixing angle solution by the KamLAND data. 

Yet, there are still uncertainties which arise due to the two remaining unknowns, the 
neutrino mass hierarchy and the value of #13. The former leads to the freedom of the H level 
crossing in either the neutrino or the anti-neutrino channels depending upon the normal or 
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Flux model averaj 


?e energy (p e ) 


average energy (y x ) 


pinching (v e ) 


pinching {v x ) 


Garching (E 9e ) 


= 15 MeV 


(E Va ) = 18 MeV 


(3 = 5 


(3 = 4, 


Fermi-Dirac Tp e = 


3.71 MeV 


T„ x = 5.13 MeV 


r) = 3.12 


r] = 1.70 



TABLE I: Reference input values of the Garching and the Fermi-Dirac flux parameters used in 
this work. The total energies and the width (-E 2 ) / (E) 2 of v e and u x are taken as follows: E\°^ = 
= 5 x 10 52 erg, p Pe = (E$ e ) / (E p f = 1.2, and p Px = {E„ w ) 2 = 1.25. The flux *„ a 

is determined by the relation £'*° t = <& Va {E Va ). The temperatures and the r\ Va of the Fermi-Dirac 
distribution are chosen to reproduce the same average energies listed in the upper row and the 
same p. 



the inverted mass hierarchies, respectively, which in principle can be used to discriminate the 
hierarchy 481.149. 
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561 ] . The latter allows the possibilities of the adiabatic 
or the non-adiabatic H level crossing depending upon sin 2 #13 > 10~ 3 or sin 2 #13 < 10~ 5 . 



III. ASSUMPTIONS 
A. SN flux model 

We employ a SN model-dependent fitting method to reconstruct SN neutrino fluxes at 
the core. We trust global features of the SN neutrino fluxes at the core obtained by detailed 
simulations and parametrize them by simple functions. We adopt two choices, namely, the 
traditional Fermi-Dirac form Eq. (CQ) and the Garching parametrization in Eq. (jSJ). Both 



parametrizations contain three parameters, (E Ua ) , (3 Ua ) and ($„ Q , T h 



20 



21 



for the 



56 



Garching and the Fermi-Dirac distributions, respectively. In previous works 
the pinching parameters were always fixed and thus assumed to be known in advance. The 
inclusion of (3 Va (r] Va ) parameter as a fit parameter constitutes therefore a new feature of 
the present analysis. We feel that this inclusion is essential because the pinching parameter 
represents the effect of departure from local thermal equilibrium and would reflect the envi- 
ronmental condition of matter around the neutrinosphere. Hence, it should be determined 
by observations. The values of initial parameters used to generate SN neutrino "data" are 
summarized in Table [H We assume that SN neutrino spectra are time-independent during 
the burst. Although this may be too idealized, it is a conservative assumption in the sense 
that relaxing this it would further complicate the task of resolving the degeneracy. 
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B. Detector and analysis method 



In this paper, we focus on water Cherenkov detectors as they are most likely the ones 
to run sufficiently long so as to watch galactic supernovae over long enough time scales. 



In particular we consider Hyper-Kamiokande 24|, and take the fiducial mass of 720 kton 



assuming that the whole inner volume is available for SN neutrinos. For other similar 



projects see e.g. Refs. |25j, |26|] 



In this kind of detectors several reactions contribute to the SN neutrino signal: inverse 
beta decay, elastic scattering off electrons, and CC and NC with oxygen. Nevertheless, 
inverse beta decay is indeed the dominant reaction, which would yield ~ 2 x 10 5 events in 



Hyper-Kamiokande for a SN at 10 kpc from the Earth. For the cross section see Ref. 581.159]. 
The expected event numbers from v e and v x elastic scattering and v e absorption by Oxygen 
are more than an order of magnitude below [3] . 

In this paper, thus, we restrict our analysis into a unique observable, the positron energy 
spectra produced by the v e absorption reaction on protons, v e + p — > e + + n. For this 
reaction only six parameters are relevant: E^, (E Pe ), p Pe = (-£^ e ) / (E Pe ) 2 , E^, (E V J), 
and p Vx = (E 2 J) / (E Ux ) 2 . We assume the SN at a distance of 10 kpc from the Earth and 
generate the SN neutrino flux data by assuming the seed parameters listed in Table [B The 
data are subsequently fitted using the same parametrization assumed to generate it, except 
for results shown in Fig. [9l Motivated by our current understanding of the composition of 
the proto- neutron star we also assume that the ratio 3 r = (E Ux ) / (E Pe ) > 1. We perform 
a standard x 2 analysis including only the statistical error. The energy bins used are chosen 
as: 70 bins for 5MeV < E e < 40MeV, 10 bins for 40MeV < E e < 50MeV, 3 bins for 
50MeV < E e < 56MeV, 1 bin for 56MeV < E e < 60MeV, and the highest energy bin 
for 60MeV < E e < lOOMeV (total 85 bins), where E e stands for the measured energy of 
the positrons. The energy resolution is taken into account with a Gaussian function with 



width a TCS = 0.47^/.E e /MeV MeV, following Ref. 60]. The results are represented in terms 
of E%, (E Va ) and (E% a ) / (E U J 2 . The 2 (3) a CL allowed regions are determined by the 
condition, 

Ax 2 = X 2 -xLn<6.18 (11.83), (6) 



3 To extract the real degeneracy we want to eliminate a trivial degeneracy which inevitably comes in into 
such fitting procedure. In view of Eq. ([7]) it is clear that for any solution for the six parameters labeled as 
(a), there is another solution labeled as (b) in which (3^, (EpJ^ are replaced by , and (E^ m )^ a \ and 
vice versa, together with = t an 2 ^ 12 $^ ; $W = co t 2 12 $£^. In these expressions, we have used the 
approximation S13 -C 1 for simplicity. We have explicitly verified that the exchange degeneracy solution 
can be removed by imposing the former condition. 
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for 2 degrees of freedom. 



C. Neutrino mass hierarchy and validity of our ansatz of SN neutrino flux spectra 



In this paper we assume that SN neutrino spectra at the surface of the progenitor star 
are given by the flavor transformed ones of the initial spectra of either the Garching or the 
Fermi-Dirac type. 

In order to illustrate the possibility of inferring the original neutrino spectra we consider 
the case with normal mass hierarchy. Ignoring the Earth matter effect, the relationships 
between the v P flux at the core of SN and the one at terrestrial detectors is as follows: 



F _ _ r 2 2 F , M _ 2 2 \ t?0 2 rpO , 2 r^O 



(7) 



where the expressions s\ 2 , c 2 2 , sf 3 , and c 2 3 stand for sin 2 #12, cos 2 #12, sin 2 #13, and cos 2 6*13, 
respectively. The observed v e spectra is then a superposition of the original v e and u x spectra. 
The coefficients of the composition depend basically on the value of #12, which we will assume 



to be sin 2 ^12 = 0.3 la, LlTl] - Therefore the goal of the analysis is to disentangle the two 



different components present in the observed v e spectrum. 

The validity of our analysis would be affected if any mechanisms are operational inside or 
outside SN core that invalidate the above ansatz. These include basically three possibilities. 
The first one is a possible unknown time dependence of the parameters characterizing the 
initial neutrino fluxes. In this case one can always carry out the analysis for considering 
the positron spectra at different time bins. The main consequence of this effect will be an 
increase of the statistical errors but the main features would remain. 

The second aspect that could affect our analysis is the presence of the shock waves 



propagating within the supernova. In several works [61|, |62j, |63|, [64| it has been discussed 
how the presence of these discontinuities could introduce an energy and time modulation in 
the survival probabilities. However, the v e flux is only affected in the case of inverted mass 
hierarchy and "large" 6 , 13 , sin 2 #13 10 -4 . The formation of shock waves could additionally 
lead to turbulent density fluctuations, causing a neutrino flavor depolarization as discussed 
in j^] and (66 ] where the former considered the 5-correlated density fluctuation whereas 
the latter considered the Kolmogorov-type turbulence implied by realistic SN simulations. 
The characteristic signal due to the shock wave, which was originally considered visible 
ignoring the effects of turbulence and/ or density fluctuations, tend to be washed out by 
such turbulent density fluctuations [65|, [66] . See also [67], |68| for effects related to the shock 
waves. We note, however, that as long as the v e flux (at Earth) is concerned, these effects 
would be significant only in the case of inverted mass hierarchy and not so small value of 
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13 , sin 2 # 13 >0(1(T 3 ). 

The third case where our analysis would not apply would be in the presence of deco- 
herence. This can be due to the multi-angle nature of the neutrino self-interaction, and 
can be present for both normal and inverted mass hierarchies. Nevertheless, it is expected 
to significantly affect the neutrino propagation only for very similar y e and v e spectra 40], 
which does not seem to be the case according to the SN simulations 32|. Therefore we will 
neglect it. 

Last, but not least, the possible existence of non-standard neutrino interactions, apart 
from affecting neutrino propagation through the SN envelope, could induce resonant conver- 
sions in the most deleptonised inner layers 69| with possibly dramatic effects in a Megaton 
water Cherenkov detector 



70| | . This possibility will also not be considered here. 



Before closing this section let us mention that a similar analysis could be performed in 
the case of inverted mass hierarchy. However, in this case one should take into account the 
possible modulations induced by the shock wave passage, as well as by collective effects such 
as the spectral split in j^J or the effect of the second-order difference between the 
and v T refractive index |43| . This comparative study lies however beyond the scope of this 
paper. 



IV. DEGENERACY IN FITTED PARAMETERS 
A. Continuous degeneracy in the fit parameters 

Let us start with the analysis with the Garching parametrization of the SN v fluxes 
assuming the values for the initial v spectra given in Tab. [J To understand the effect of 
varying pinching parameters let us first assume that the f3 Va parameters are known. In Fig. H] 
we represent by the dashed (solid) blue ellipses the regions allowed at 2a (3a) CL in the 
space spanned by (E Pe ) — {E Vx ) (left panel), {E Ue ) — Elf (middle panel), and (E Vx ) — E l °* 
(right panel). In each panel, the best fit point is also indicated by a star, which of course 
reproduces the input value. It can be seen how we can determine (E Pe ), (E Ux ), E^, and E£ ot 
with an accuracy of roughly 2%, 4%, 15%, and 30%, respectively, at 3a CL. These results are 



in good agreement with our previous work 20J apart from small differences due to different 



assumptions on the initial spectra and the detector (now with slightly smaller volume). It is 
remarkable that the non-vanishing mixing angle 9 12 allows us to obtain information about 
v x flux parameters event though only z/ e 's are directly detected. This is a direct consequence 
of Eq. ©. 

However, once we allow the pinching parameters (3 Va to vary freely, the accuracy in the 
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FIG. 4: Comparison between the 2a CL and 3a CL determination of the astrophysical parameters 
for the cases of fixed and free pinching parameters, displayed in terms of (Ep e ) — (E Ux ) (left panel), 
{Ep e )—Ep°^ (middle panel), and (E Ux )—E v ° x (right panel). The solid and the dashed ellipses indicate 
the cases in which (3p e and f3 Vx are fixed to (3p e = 5 and (3 Ux = 4. The star represents the best-fit 
point, which, in this case, coincides with the input values. The constraint r = (E Ux ) / (Ep e ) > 1 has 
been imposed. The corresponding at 2a CL and 3a CL determinations for the case of free pinching 
parameters are denoted by the (dark/red) and (light/yellow) regions. See the text and Tab. U for 
more details. 



determination of the flux parameters is significantly reduced. This is represented in Fig. [4] 
with shaded areas, yellow (light) and red (dark) corresponding to 2 and 3a CL respectively. 
The striking consequence is the emergence of a continuous parameter degeneracy 4 : there 
is a continuum of allowed fit solutions. For definiteness, in what follows we make a very 
conservative assumption, (Ep x ) < 32 MeV. We find that this degeneracy affects v e and v x 
flux parameter determination in a different way. For v e the sensitivities to (E De ) and are 
reduced to 4% and 50%, respectively, at 3a CL. As can be seen in Fig. [4] for the case of u x , 
however, the effect is much more drastic. The region consisting of degenerate solutions forms 
a quasi one- dimensional strip extending mainly in the direction of v x SN flux parameters. 
This is a direct consequence of the different weights with which the original v a fluxes enter 
in the observed v e flux: 70% from v e and 30% from u x , see Eq. (JTj). 

In order to shed more light on the nature of the degeneracy we present in Fig. [5] the 
allowed regions displayed in terms of various other combinations of the flux parameters. 
^From the last panel of Fig. H] and the top left panel of Fig. [5] it can be inferred that the 
key of the continuous degeneracy lies on the possibility of increasing the pinching of F® , 



Note that this degeneracy is quite different in nature from the (discrete) degeneracy one encounters in 
the determination of lepton mixing parameters in neutrino oscillation experiments [7ll. |72|. |73|. 
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1,1 1.15 1.2 1.25 1.3 1-15 1.2 1.25 1.3 1.15 1.2 1.25 1.3 

<e 2 (k)>/<e(k)> 2 <E 2 {V,)>/<E{v,)> 2 <E 2 (u.)>/<E(F,)> 2 

FIG. 5: Best-fit point (star) and 2a CL and 3d CL projections over the displayed parameters for 
free pinching parameters in the Garching parametrization. 

i.e. larger (5 Ux and smaller p Ux . By requiring a stronger pinching of the v x flux, together 
with a larger value of (E Vx ) and a simultaneous reduction of i?* ot it is perfectly possible to 
mimic the behavior of observables defined with input values of F® . In the bottom right 
panel one can see how the required variation of p Vx does not necessarily imply a significant 
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modification of the shape of the v e flux. This feature can be seen for instance in the left 
panel of Fig. [TJ In that example the fake spectra can be made surprisingly similar to the 
true one roughly by keeping the same original F® parameters and changing F° following 
the previous recipe: a reduction of E^ together with an increase of (5 Vx and (E Ux ). 

Another salient feature that can be observed in most of the panels of Fig. [5] is an island 
structure, i.e. the presence of a region separated from the main allowed region containing 
the best fit point. To understand the origin of such structure one has to realize that this 
region arises for larger values pp e and (E Pe ) than the initial ones, and smaller p Vx and (E Vx ). 
This is related to the trivial degeneracy mentioned before, the extra solution obtained by 
interchanging the v e and v x spectra. When we do not impose the condition (E Ux ) > (Ep e ), 
the allowed region has a "cross shape" because of the trivial degeneracy, namely the region 
with swapped parameters between v e and v x is allowed. It appears to us that the island 
which is left over is a remnant of the swapped parameter region. 

B. Robustness of the degeneracy 

It is reasonable to ask whether the degeneracy is an artifact of the particular parametriza- 
tion we employ, or an accidental consequence of the particular choice of the initial param- 
eters. In order to answer the first question we performed the same exercise by using the 
pinched Fermi-Dirac parametrization, Eq. (0Q), both in the preparation of data, using the 
input values given in Tab. [H as well as in fitting them. As shown in Fig. [6] the shapes and 
the sizes of the allowed regions are very similar to those given in Figs. H] and [5] (for brevity, 
we have shown our results only for 6 representative combinations out of 15, because the 
others are similar to previous plots). The continuous parameter degeneracy is also present 
for the same combination of v x parameters: strong pinching, large mean energies and small 
integrated luminosity. The main difference from the case with the Garching parametrization 
is the absence of the island structure. 

In order to shed some light on the dependence of the results on the values of the initial 
neutrino fluxes we have redone the fit in the case of Garching data fitted by Garching 
distribution for (E Ux ) = 16.5 and 20 MeV, which are presented in Figs. [7] and El respectively. 
One can see the shape of the allowed regions changes but the presence of the continuous 
degeneracy persists with a similar range of v x parameters. On the other hand, the size of 
the island region changes depending on the initial parameters considered. The more similar 
the initial v e and u x spectra, the bigger the island. This is due to the fact that it is easier 
to interchange the role of the two flavors. 

For completeness we have also varied the number of energy bins used in the statistical 
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1.1 1.15 1.2 1.25 1.3 1.1 1.15 1.2 1.25 1.3 1.16 1.17 1.18 1.19 1.2 1.21 1.22 

<E 2 (k)>/<Uk)> 2 <E. 2 (k)>/<E(k)> 2 <E 2 (i7 e )>/<E(i7,)> 2 



FIG. 6: Best-fit points (star) and 2a CL and 3a CL projections over the displayed parameters for 
free pinching parameters in the Fermi-Dirac parametrization. Here six representative combinations 
of fitting parameters were chosen, corresponding to the 3 panels in Fig. |4] and top 3 panels in Fig.0 

analysis, but found no significant change in the result. We conclude then that the presence 
of the continuous degeneracy is a robust feature of the reconstruction analysis of supernova 
parameters. 

C. Effect of uncertainties in the SN neutrino flux spectra 

So far we have assumed in our analysis that we know the functional form of the supernova 
neutrino flux spectra prepared by the exploding star. Of course, this is not the case. What 
would be the effect of our ignorance of primordial SN neutrino flux spectra parametrization 
on the analysis? In order to gain some insight on this issue in this subsection we attempt a 
new procedure: to generate the data with a Fermi-Dirac distribution and to fit it assuming 
the Garching parametrization, and vice versa. This way we try to check not only whether 
it is possible to determine the flux parameters, but also its functional form. 
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14 14.214.414.614.8 15 15.215.4 

< E(F„) > (MeV) 



1.15 1.2 1.25 1.3 
<E\u x )>/<E(p,)> 2 



FIG. 7: Best-fit point (star) and allowed 2a CL and 3a CL projections over the displayed parame- 
ters for free pinching parameters in the Garching parametrization. Here we use as input the values 
given in Tab. [U except for {E Ux ), which is 16.5 MeV. Two fitting parameter combinations were 
selected: the left panel corresponds to the top left panel of Fig. [6] whereas the right one corresponds 
to the bottom left panel of the same figure. 




.6 14.8 15 15.2 15.4 15.6 15. 

< E(P.) > (MeV) 



1.15 1.2 1.25 1.3 
<E 2 (v«)>/<E(0> 2 



FIG. 8: Same as in Fig. [7] but for (E Ux ) = 20 MeV. 



The results of such a fit are presented in Fig. In the left panels we show the allowed 2a 
and 3a CL contours assuming data generated with a Fermi-Dirac distribution, with input 
values given in Tab. [J and fitted with the Garching parametrization. First of all, one realizes 
that the shape of the contours are very similar to those presented in the previous sections, 
see e.g. Figs. H] and [51 In particular, one finds the same one-dimensional correlations between 
the different parameters characterizing F® . It is curious to note that the size of the allowed 
regions is smaller than in the previous figures. The second feature to point out is that the 
best-fit point does not coincide with the initial values assumed to generate the data. The 
X 2 value obtained in the best-fit point is xlets = 5.29 5 . 

5 Notice, however, that these numbers have mild dependence on how fine is the mesh by which the parameter 
space is covered. 
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1.1 1.15 1.2 1.25 1.3 1.1 1.15 1.2 1.25 1.3 

<E 2 (i/ x )>/<E(0> 2 <E 2 (l/ x )>/<E(l/,)> 2 




< E(>„) > (MeV) < E(z/„) > (MeV) 



FIG. 9: Left panels: Allowed 2a CL and 3a CL contours assuming data generated with a Fermi- 
Dirac distribution and fitted with the Garching parametrization. Right panels: Allowed 2a CL 
and 3a CL contours assuming data generated with the Garching distribution and fitted using the 
Fermi-Dirac parametrization. The initial values are taken from Tab. U and are represented by a 
filled circle. The best fit points are shown with a star and have a value of Xbest = ^-29 and 3.34 in 
the left and right panels, respectively. 

In the right panels we show the same contours fixing the data with the Garching 
parametrization and carrying out the fit assuming a Fermi-Dirac distribution. One can 
still recognize roughly the same tendencies as in the previous case. However the contours 
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are much bigger. As far as the best fit point is concerned, it has a value of XtLts = 3.34, 
and as in the previous case it also does not coincide with the input values. These broad 
features indicate that Fermi-Dirac parametrization has somewhat higher flexibility than the 
Garching one, at least for the input values we used in this paper. 

We note here that the results presented in Fig. [9] not only demonstrate the robustness of 
the degeneracy but also serve for providing some clues to answering the question of whether 
the SN neutrino observation can discriminate between different energy distribution of SN 
flux. The values of xLst f° r both cases (input given by Garching parametrization fitted by 
Fermi Dirac one, and vice versa) are small compared to the number data points (bins), 85 
bins. This implies that both distributions can fit equally well irrespective of the assumed 
input. Within our current procedure it seems rather difficult, by fitting only the observed 
data, to say which parametrization would describe better the SN spectra. We note that 
such a difficulty arises because the observed spectra are given by the superposition of the 
2 spectra which is the origin of the degeneracy we found in this work. If the spectra were 
described only by the single spectrum by either Fermi-Dirac or Garching parametrization 
(hence no degeneracy) it would be possible to distinguish the functional form thanks to the 
huge statistics. 

D. Analytic understanding of correlations 

Examining previous figures, e. g. Fig. HI Fig. [5] and Fig. El one sees that they exhibit 
some intriguing correlations between the flux parameters. Let us now try to understand 
these features. We have seen in Fig. [1] that the degeneracy is present at the level of the 
fluxes arriving at the Earth, already before the detection process enters into the game. In 
particular, since our detection channel is the inverse beta decay, all the features observed 
in the fit should be encoded in the v e flux at the Earth. For this reason let us consider the 
following quantities 

Eg = J EF„.dE t (E oh ) = J j E ^f^ (E 2 oh ) = J f F F p f E E (8) 

These functionals are determined in terms of the energy spectrum of F Pe at the Earth. 

Of the six fit parameters we have, one notices from the previous Figures that the range of 
variation of two of them is rather narrow, namely {E Pe ) and pp e . Hence in order to understand 
the nature of the correlations we need only the 4 remaining "effective" parameters. Using 
Eq. ([HD and requiring these quantities to be equal to the observed (input) values, we see 
that these 4 parameters are correlated, and one can choose to express the dependence of 
(E Ux ), and p Vx , in terms of one, e.g. E x °\ 
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Taking into account the relation in Eq. (J7|), is trivially obtained by our assumption 
on the initial values of the total energies = E^ = 5 x 10 52 erg, see Tab. [H as 

E£ = <M? + 42K o : = 5 x 10 52 erg , (9) 

where c 2 2 and s\ 2 stand for cos 2 9\2 and sin 2 # 12 , respectively. If we now fix E^ to this value 
and let E^ and E^ freely vary we get, after substituting the value of 6*12 the following 
trivial relation 

Elf = - 0A3E™ , (10) 

which explains the anti-correlation found between El° l and i?* ot in the right panel in the 
second row of Fig. [5j 

The approximate quasi one-dimensional correlation is displayed in Fig. [TO] In the figure, 
we plot a = a(E^), where a represents El° l , {E Pe ), {E Ux ), p Pe , and p Ux . It is defined as the 
value of a that minimizes x 2 { a i -^vf ) f° r a given E l °\ This plot illustrates how the different 
parameters adjust themselves when i?* ot is varied so that the new set of parameters fits 
equally well. Together with a(E^ >t ) we show the analytical expression obtained by requiring 
the quantities defined in Eq. (JSj) to be equal to the observed (input) values. 

In the top panel of Fig. [10] we can see how the expected anti-correlation for a = E^ 
(see the right panel in the second row in Fig. [5]) completely agrees with the analytical 
expression [101 

Let us now consider the mean energy of the v e arriving at the Earth. According to the 
values assumed for the input parameters describing the flux one expects, 

ptot 

< ^ >= w^)+ w ^) =15 - 79Mev - 

We can now proceed as before and rewrite (E Ux ) in terms of E^ taking into account Eq. (Q: 

{EuJ = % 2 E%-E%[l- <^ e >/(£ob>) (jEpe) = ™ E£- 8.34x10^ erg MeV " (12) 

where in the last step we have assumed that (E Pe ) and (E D b) are constant and equal to 15 
MeV and 15.79 MeV, respectively. Under this assumption we expect that (E Ux ) decreases 
for large values of E^. This behavior can be clearly seen in the right panel in Fig. [4] as well 
as in the middle panel of Fig. [10] (see the red solid curve). Moreover, from the same panel 
one sees that the assumption of constant (E Pe ) (blue curve) is justified when compared to 
the variation of (E Ux ) (red curve). However, this assumption is not perfect and therefore the 
observed variation of (E Ux ) and the analytical prediction from Eq. [12] (black curve) do not 
coincide. 
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FIG. 10: Various quantities are plotted as a function of , in blue and red for v e and u x , 
respectively. The total energy E}f is shown in the top panel, the average energies {E) in the 
middle panel, and the width parameters (E 2 ) / (E) 2 in the bottom one. The analytical values 
predicted in Eqs. (fTU| PT2| [Til) are also plotted in black. See text for details. 



Now we analyze the second moment of the flux, 

{E ° h) ~ c 2 2 E^/(E Pe ) + s 2 2 E%/(E Ux ) " 305 - 53 MeV - ( 13 ) 

Taking into account Eqs. (GO), ( TTTj) and ( jl~3l) we can express p Vx in terms of (E Ux ): 

„ (E 2 J f / (E Pe ) \ (E, e )/(E oh ) - (Ej)/(E 2 h ) (Ej) \ 

VUX (Eo a ){E Vx )\\ (E Vx )j l-(E, e )/(E oh ) + (E 2 J\- 1 ] 

Analogously to the previous case we can require (E 2 h ) and p„ e to be constant. Then, from 
this expression one can see how as (E Vx ) increases the v x spectrum has to become more 
pinched in order to keep the same shape. Taking into account that (E Ux ) and E^fJ 1 are anti- 
correlated one expects a positive correlation between p Vx and i£* ot . This is exactly what 
is seen in Figs. |5] (left panel in the second row) and in the bottom panel of Fig. [10] (red 
curve). In the latter figure we see how this argument is mainly qualitative, since p Pe is not 
constant over the whole range of the E^ considered. Actually, when comparing this result 
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with the analytical condition (fT4j) one realizes that this qualitative argument only applies 
for a definite range of . 

V. RESOLVING THE DEGENERACY? 

Once the presence and robustness of the degeneracy have been established, the next step 
is to point out strategies to solve the problem. Here we comment on three possible directions 
that could be followed. 

First from the theory side we note that a better understanding of the neutrino flux 
formation would narrow down the range of values that the flux parameters could take. In 
particular we have seen in Fig. [5] how our knowledge of the pinching parameter can play an 
essential role in reducing the uncertainties on the traditionally considered parameters (E Va ) 
and El ot . 

From an experimental point of view one may correctly argue that the degeneracy we 
have pointed out is an artifact of our treatment, which uses only the v e absorption reaction. 
The obvious way of resolving such a degeneracy would be to include other reactions into 
the analysis, such as neutrino elastic scattering off electrons, or CC and NC reactions with 
oxygen. Indeed, the quasi one- dimensional nature of the degeneracy we found would suggest 
that it would be lifted by adding an extra high statistics observable, independent of inverse 
beta decay. 

However, the situation in our case is quite different. To understand why, one must recall 
the following crucial fact. In order to include the previous reactions one must also include 
v e into the analysis. As a result instead of the six supernova neutrino flux parameters used 
in our fitting procedure we would need nine. Since the number of events due to the above 
reactions is more than one order of magnitude less than those coming from the v e absorption 
reaction, it is not obvious at all that, by adding these other reactions one will indeed resolve 
the degeneracy we have uncovered. 

If the degeneracy cannot be resolved simply by using water Cherenkov detectors one must 
think of possible alternative detectors to combine with. For example the situation might 
improve either by: (1) adding detectors with better sensitivity on v x parameters, or (2) 
adding detectors sensitive to v e . Good candidates in the current and near future experiments 
which can contribute towards these goals would be Borexino (74)], KamLAND (in its future 
solar neutrino observation phase) 75] and HALO project The former 2 detectors will 



be able to observe proton recoil in vp elastic scattering in a liquid scintillator 77J. We 
stress that this would provide a unique opportunity to obtain the spectrum information of 
NC reactions which should be very important in resolving the degeneracy. The latter, on the 
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other hand, would have capabilities to detect v x and v e through the NC and the CC reactions 
on lead. With regard to the possibility (2) above it would be interesting to consider a high 
statistics measurement of v e events in a liquid argon neutrino detector through the charged 



current absorption of u e by 40 Ar 57J. As can be inferred from Eq. (JTj), the dependence of 
the v e flux at the Earth on the original F® and F9. is different from that of antineutrinos. 
In particular for large #13 the following condition holds 

*u = 4,*£ + (is) 

Therefore, one could hope that this fact could break the degeneracy expected in water 
Cherenkov detectors. Nevertheless the possibility of distinguishing between the u e and v e 
signals is in any case non-trivial. Therefore, a realistic nine-dimensional parameter analysis 
taking into account all neutrino flavors is outside the scope of this paper but should be taken 
up as a next step. 



VI. CONCLUDING REMARKS 

In this paper, we have reexamined the possibility of reconstructing the initial fluxes of 
the neutrinos emitted in a future galactic core-collapse supernova by using a Megaton-sized 
water Cherenkov detector. 

The three parameters that are usually considered to characterize the non-thermal su- 
pernova neutrino flux are the average and total energies of each species, and the so-called 
pinching parameter. The latter is connected to the second moment of the distribution func- 
tion and modulates its shape. Due to the current uncertainty on its precise value, we have 
included the pinching parameter as a fit parameter, in contrast to previous works where 
only the average and total energies were considered. 

In order to illustrate our results we have considered the following scenario. First of all, 
from all reactions giving rise to the neutrino signal in a water Cherenkov detector we have 
concentrated on the inverse beta decay, v e +p — > e + + n, which is by far the most important 
one. Therefore, our setting is sensitive only to the anti-neutrino fluxes at terrestrial detectors. 

As far as the neutrino properties are concerned, we have focused on the case of normal 
mass hierarchy neutrinos, neglecting also any non-standard neutrino interactions which could 
affect neutrino propagation in an important way. Such well-defined scenario is much less 
affected by uncertainties due to neutrino self-interaction in the inner layers or the distortion 
of the matter density due to the shock wave passage. 

On the other hand the v e flux arriving at the Earth is in this case a strong mixture of the 
initial v e and v x . Therefore the detected v e flux would provide us with information not only 
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about v e but also v x . Under these assumptions we performed a \ 2 analysis on artificially 
generated data. 

We found that the inclusion of the pinching parameter in the fit analysis has a drastic 
consequence: the appearance of a continuous degeneracy in the determination of the v x flux 
parameters. This degeneracy is quite robust, as it persists irrespective of the parametrization 
taken for the neutrino distributions. This makes very difficult a complete determination of 
the different parameters characterizing the v x flux using the inverse beta decay events, even 
in the case of a Megaton water Cherenkov detector. 

The solution to this degeneracy must come from a better understanding of the neutrino 
spectra formation, as well as an optimization of the information provided by the comple- 
mentary neutrino reactions involved not only in water Cherenkov detectors, but also in 
alternative detector techniques. 

Note added: After completion of our work we became aware of the paper 78| . which was 
in fact triggered by previous private communications with us. Their analysis includes other 
reactions available in water Cherenkov detectors. However, their results indicate that the 
addition of these channels is not enough to fully resolve the degeneracy pointed out here. 
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